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Abstract 

We study the distribution of the complex temperature zeros for the 
partition function of the Ising model on a Sierpinski gasket using an exact 
recursive relation. Although the zeros arrange on a curve pinching the 
real axis at T = in the thermodynamic limit, their density vanishes 
asymptotically along the curve approaching the origin. This phenomenon 
explains the coincidence of the low temperature regime on the Sierpinski 
gasket and on the linear chain. 



1 Introduction 

The understanding of phase transitions on non-crystalline structures has been 
recently improved by exact results connecting general geometrical features of 
networks to the existence of spontaneous symmetry breaking jl|, 0, 0, |J. The 
situation is more complex when dealing with the critical behaviour of model sys- 
tems. For continuos symmetry models the singularity of the free energy, which 
determines the critical behavior, appears to be related to the infrared spectrum 
of the Laplacian operator on the network ^| , while this is not the case for dis- 
crete symmetries (e.g. the Ising model). There, all known results suggest that 
the link between critical behaviour and geometry should involve some other 
topological features 0, 0]. An interesting result concerns the Sierpinski Gasket, 
a typical and widely studied fractal, where the Ising model is exactly solved. 
On this structure, although continuous symmetry models exhibit a power law 
behaviour for T — ► 0, the Ising model has an exponential low temperature be- 
haviour which coincides with that found on the linear chain pi 0] . To analyze 
the critical regime from an ab initio point of view, an interesting picture is pro- 
vided by the study of the singularities of thermodynamic potentials. In 1952 Lee 
and Yang in two famous papers [0 ^| first proposed their fundamental approach 
to phase transitions, consisting in studying the zeros of the partition function 
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of a statistical system, considered as a function of a complex parameter. The 
partition function on a finite volume is a polynomial in complex activity or fu- 
gacity, so that the complete knowledge of the zeros distribution is equivalent 
to the knowledge of the partition function itself and all thermodynamic quan- 
tities can be obtained from it. On a finite volume there are no real zeros, the 
coefficient being all real and positive. However in the thermodynamic limit the 
zeros can pinch the real axis (the region of physical interest) producing a singu- 
larity in the free energy (or grand-canonical potential) [[)) . The pinching points 
are phase transitions points on the parameter axis and the zeros distribution in 
their neighbourhood can be connected with the critical properties of the system 
& 

Unfortunately, a complete knowledge of the zeros is very difficult to obtain 
except for a few exactly solvable cases. General theorems hold for the zeros 
distribution in the complex magnetic field plane in a class of ferromagnetic 
lattice systems, including the Ising model [§> f f . On the other hand, very little 



is known rigorously about the behaviour of the zeros of the partition function 
in the complex temperature plane, the so called Fisher zeros p2[ . In general, it 
is not clear if Fisher zeros arrange on smooth curves even if this is the case in 
some exactly solvable models. For the Ising model on regular two dimensional 
lattices , the Fisher zeros arrange on curves that cross the positive real 

axis at the transition point. In the one dimensional case only two zeros (with 
infinite multiplicity) are found and these have a nonzero imaginary part, so that 
there is no singular point for the free energy. 

For statistical models defined on non periodic discrete structures, Fisher 
zeros show some peculiar features, making the analysis of their density and 
location extremely subtle. In particular, on some hierarchical lattices (i. e. 
q-potts model on diamond hierarchical lattices [|l4|) the zeros have been show 
to form a fractal set (Julia sets). In this case, while the general approach 
for identifying the singularity points and the critical behaviour still holds, the 
widely used arguments concerning scaling of singularities and zeros density with 
the volume must be handled carefully, as will be shown in the following. 

In this paper we will study the Ising model on the Sierpinski gasket, obtaining 
a recursive relation for the partition function, from which the zeros of the n-th 
stage gasket can be obtained from those of the (n — l)th. The distribution we 
obtain is fractal and pinches the real axis at T = 0, so that a singular point 
with a power law critical behaviour could be expected. However, since the zeros 
density is found out to vanish exponentially in the neighbourhood of T = 0, 
these zeros don't produce any singularity of the free energy: although the zeros 
pinch the real axis the 'critical behaviour' is the same as the one-dimensional 
case. 
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Figure 1: First iterations of gasket's construction 

2 Ising model on the Sierpinski gasket 

The Sierpinski gasket is a fractal graph which can be built recursively with the 
following procedure: the initial stage (Go) is a triangle (3 sites with 3 edges) 
and the nth stage (Gn) is obtained joining 3 Gn-i at their external corners, to 
form a bigger triangle (fig. [j]). In this way Gn has |(3™ — 1) sites, 3™ +1 edges 
and its side contains 2™ edges. 

The gasket is obtained as the limit for n — > 00 of this procedure. 

The Ising model on the gasket is defined associating the spin variable <7j = ±1 
to every site i of the graph, and considering a nearest-neighbours interaction 
between points joined by an edge (link). The Hamiltonian is therefore: 

E = - J E a ^ C 1 ) 

<i,j> 

where the sum runs over the couples of sites joined by a link and J is a positive 
constant (ferromagnetic coupling). 

3 Recursive relation for the partition function 

For Go the partition function 

Z=Y. ( 2 ) 

{^} 

can be seen as a sum of the elements of the rank 3 tensor Mq 

Z Q = M) CT1<T2CT3 ( 3 ) 
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Figure 2: Labeling of sites used for Qq and Q\ 



where 



M w°> = exp{-[3E(a u * 2 ,<T 3 )} (4) 

Mq 1 ' 72 ' 73 can take only 2 values because there are only 2 classes of spin config- 
urations with different energy: 

71 fooo _ 3PJ _ ,,3 

M - e — V fe-\ 
M™ ( ~ a) = M° { -° ] ° M^ a)aa = e-? J = y- 1 [ ' 

where y — e@ J . In terms of y the partition function is: 

Z a = 2y 3 + Gy- 1 (6) 

For Q\ the partition function can be expressed in the same way separating 
the sum over the states of the inner sites using the tensor Mi whose indices 
correspond to the spins on the external vertices (fig. ||): 



M vi°2<T3 = J2 exp[-PE(<n)] = J2 M^^M^^M^"" 3 

(T a ,<7b ,(T c = il £7 a ,<7f> .£7 C = ±1 

(7) 

Now Mi has the same structure as Mq since the possible values are: 

M ™(-°) = M ^- a)a M[- a) ' 7a = 3y- 3 + 4y + y 5 (8) 

One can obtain M\ from Mo simply by a transformation mapping y 3 in 4y -3 + 
3j/ + y 9 and y~ x in 3y~ 3 +4y + y 5 . This is done by the substitution 



y 4 + 3 



followed by the multiplication by 



c(y)-^^[( ; V 4 + 3) 3 ( 2 / 8 - 2 / 4 + 4)]' (10) 
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The transformation also gives the new partition function: 

Z 1 {y) = Z Q {f{y))c{y) (11) 

Following the same argument one can obtain for the (n + l)th stage of the 
gasket Gn+V 

z n+1 (y) = z n (f(y))-Hy)f (12) 

Using this recursion relation we get: 

Z n (y) = 4^(J/ 4 ) (13) 

yd 

where P n (t) is a polynomial in t of degree 3™ in which the coefficient of t 3 is 
1; for n — one has Po(t) = t + 3 while the general case n > can be proven 
by induction. 

4 Zeros of the partition function 

Introducing the variable x — y 4 the transformation (0), ( |To| ) is given by: 



(14) 



\ c(x) = {x + l)x-i [(x + 3) 3 {x 2 - x + 4)] ^ 
Denoting by x l n the zeros of P n (x) the partition function reads 

Z n = Jqi fl(^ - O (15) 
i— 1 

and using the recurrence one finds 

+i 3 " +1 +i 3 " 

2x" a V Y[ (x - < +1 ) = 2x-^ Y[{[(x 2 -x + 4)- xi(x + 3)}(x + 1)} 

i=l i=l 

(16) 

This equation shows that for every root x l n of Z n , Z n+ \ has the root x = — 1 
and the 2 solutions of 

4 = /«+i) (17) 

namely the preimages of x\ L by the transformation / . 

Starting from Xq — — 3 one obtains all the zeros of the partition function 
for the n-th stage gasket as shown in table |l|, where h{x) denotes the set of the 
preimages of x (and h k [x) is a set of 2 k zeros) . From this table one can see that 
the preimages of —1 by the jth iterate of / appear at the (j — l)th stage and are 
zeros of all the following stages, while the preimages of —3 are 'temporary' zeros. 
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n 


zeros 





-3 


1 


-1 ft(-3) 


2 


-1 ft(-l) /i(/i(-3)) 


3 


-1 ft(-l) h(h(-l)) h(h(h{-3))) 


n 


-1 ft(-l) ... h n ~ l (-l) h n (-3) 



Table 1: Zeros of partition function in y 



It is also possible to find the multiplicity of these zeros: in fact since every root 
generates the root —1 at the next stage, this value appears 3™ -1 times among 
the zeros of n-th stage, the roots h(—l) appear as many times as —1 in the 
previous stage (their multiplicity is 3™ -2 ) and in general the multiplicity of the 
roots belonging to h?{— 1) is 3 n— 

In this way one can, in principle, calculate all the zeros of the partition 
function at any stage. In practice this is possible only for small n because of 
their exponential growth. An alternative approach |l4[ ] is to start from a root 
xq (for example '—1' whose preimages are 'permanent') then choose at random 
one of its two preimages by the transformation / (denoted by xi), then choose 
one of the preimages of x\ and so on; the set of points obtained in this way is 
a representative of the set of all roots and has the advantage to contain zeros 
relative to large n. 

By plotting in the complex temperature plane the roots obtained with both 
methods it can be seen that very few of them fall in the neighborhood of the 
real axis and no information can be obtained about the critical behaviour (see 
fig. H where the zeros are plotted in the plane of the variable t = e~@ J ). 

A good technique to obtain more zeros near the real axis consists in changing 
the choice probability of the two preimages fli|| ; the two solutions of x = /(#') 
are 

x =hi(x) = (18) 



, , , , l + x + ^-15 + Ux + x 2 , . 

x = h 2 (x) = (19) 

and one can see that the repeated application of hi gives a sequence of points 
approaching the real axis. The set of roots obtained by increasing the probability 
of choosing the second preimage is not a representative set af all roots (it doesn't 
show their density not even approximately) but gives us a chance to observe their 
behaviour in the interesting area (see fig. 0). 

A plot of these roots in the plane of w = e~ 4/3J (T = corresponds to w = 0) 
with a log-log scale (fig. ||) shows that the real and imaginary part are related 
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Figure 3: 20000 zeros obtained by the random method (in the plane of t = e /9J ) 
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Figure 4: 20000 zeros obtained with probability 0.98 of choosing hi (in the t 
plane) 
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Figure 5: Log- log plot of zeros (in the variable w — e 



by a power law: the curve can intersect the real axis in the thermodynamic limit 
only at the origin. 

Analytically one can verify this power behaviour by studying the transfor- 
mation of the variable w = e _4/3J = a; -1 



fix) 



w(3w + 1) 



(20) 



_ i_ 4w 2 — w + 1 

Assuming that 3(u>) = A %t(w) b , that is 

w = £ + iA£ b (21) 
and inserting ( pl| ) in (^0|) one obtains 

+ iA£ b )) = A£ b (l + 8£ + 0(^(3,26-1}^ ( 22 ) 



and 

A (K(g(Z + iAZ b ))) b = A£ b (l + 4bZO(r in{3 ' 2b ~ 1} ))) (23) 

Choosing b = 2, one sees that the curve w = £ + iA£ 2 is 'conserved' by trans- 
formation g except for higher order terms in £. 

5 Density of zeros 



We have seen that the zeros pinch the real axis only at T = and we proceed 
by studying their density in the neighbourhood of this point to establish the 



Figure 6: Density of zeros vs. their real part in the t plane (I method) 

critical behaviour; it is important to sec whether this density (which is quite 
small, as we have seen) goes to zero or remains finite. A numerical estimate can 
be obtained by simply counting the zeros (with their multiplicity). 
This has been done in two ways: 

• considering only the zeros in the neighbourhood of the real axis and group- 
ing them with regard to their real part (one-dimensional density); 

• dividing the complex plane in equal rectangles and counting the zeros 
contained (two-dimensional density). 

In the first case one obtains plots like fig. [| which shows the density of the 
zeros of the partition function with 5(i) < .3 for gaskets from Q w to Q 2 q', in the 
second case obtains the result shown in figures |7] and || (which refers to Gis)- 
From fig. ^| one can see that going from one stage to the next the density does 
not change appreciably except for the tail towards that grows longer but is 
strongly decreasing: the density at T = appears to vanish exponentially. 

This behaviour can also be verified by an analytical estimate. First we notice 
that the zeros near T = are those obtained by the repeated application of hi. 
Indeed for |x| — > oo we have 

hi(x)->k (24) 

while 

h 2 (x) ^ x + 4 - — + 0(x- 2 ) (25) 




Figure 8: 3-dimensional view of the zeros density in the t pi 
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and, in terms of real and imaginary part, 

, . . , / , 16u \ ( 16w . 
h 2 {u + iv)tt [u + 4 +1 t) g (26) 



U 2 + V 2 J \ U 2 + V 

Applying h 2 to z = u + i v for large u one obtains h 2 (z) ~ 4m + i v. In this 
limit the density of zeros in the x plane becomes the product of two factors, one 
depending only on u and the other on v: 

d(«,w)»di(tt)da(«) (27) 

where d 2 (v) is bounded. 

The asymptotic behaviour of di(u) for u — > 00 can be obtained by noting 
that, for each set /i fc (— 1), the zeros with real part u are those obtained by a 
sequence ending with the application of h\ followed by , where n = j + c and 
c is a constant independent of u. So this fraction of zeros is | to the power ^ + c, 
that is proportional to exp(— u/U) with £7 = 4 In 2. Since the total density d\ is 
a weighted sum of the partial densities that have the same behaviour we have: 

oc exp(-u/U) (28) 

To find the density in the t plane we must now divide d by the Jacobian of 
the transformation: 

t = x~i (29) 



This Jacobian turns out to be 



±(u 2 + v*)-i (30) 
and finally 

d(t r ,U) oc d 2 {v)e- u/u (u 2 +v 2 )i ^ ^ (31) 

where t r e U are the real and imaginary part of t. 

For t — > (that is u — > +00) the density vanishes exponentially, as we could 
infer from numerical calculation, and this behaviour has the same effect as a gap 
near the real axis. Therefore the low temperature regime is not affected by the 
zeros contained in this region and one observe a situation analogous to the one 
dimensional case. This can be seen, for example, by comparing the behaviour of 
thermodynamical quantities: figure || shows that, even if the zeros distributions 
seem to be quite different, the behaviour of the specific heat for the Sierpinski 
gasket and the linear chain is essentially the same. 

6 Conclusions 

The anomalous behaviour of the density of zeros for the Ising model on the 
Sierpinski gasket is to be deeply related to the its self-similar geometry. This 



11 




Figure 9: Comparison between the zeros distribution in the t plane and 
specific heat for the linear chain (on the left) and the gasket (on the right) 
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strongly suggests a careful approach to the analysis of scaling of zeros density on 
fractals. In particular a stimulating open problem is the relation of this scaling 
with the geometry of a generic self-similar structure and with known anomalous 
dimensions. An important step in this direction would be the study of Fisher 
zeros on the more complex case of a Sierpinski carpet, where an exact solution 
is still lacking but the Ising model is expected to have a phase transition at 
finite temperature. 
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